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George  Philip  Rosser,  Jr.,  B.S.,  University  of  Georgia 
Chairman  of  Advisory  Committee:  Dr.  Phanindramohan  Das 
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Gridded  u  and  v  wind  components  derived  from  dual -Doppler  radar 
data  were  utilized  to  study  the  behavior  of  the  perturbation  pressure 
field  when  different  amounts  of  smoothing  were  applied  to  the  wind 
field.  Data-void  areas  in  the  wind  field  were  filled  by  extrapolation. 
Three-point  and  nine-point  smoothing  operators  were  employed  to  smooth 
the  wind  field  prior  to  retrieving  the  perturbation  pressure  field. 
Smoothing  the  wind  field  produced  a  nonlinear  variation  of  the  central 
pressure  value  in  a  meso-low  with  respect  to  the  smoothing  parameter. 

The  central  pressure  value  decreased  as  the  smoothing  parameter  in¬ 
creased.  Pressure  centers  less  than  0.1  kPa  (1  mb)  in  absolute  magni¬ 
tude  tended  to  amplify  and  migrate  slightly.  False  pressure  centers 
created  by  the  data-void  filling  technique  were  removed  when  a  smoothing 
parameter  greater  than  or  equal  to  0.4  was  employed.  Smoothing  para- 
meters  greater  than  or  equal  to  0.5  created  a  false  pressure  center. 
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Gridded  u  and  v  wind  components  derived  from  dual -Doppler  radar 
data  were  utilized  to  study  the  behavior  of  the  perturbation  pressure 
field  when  different  amounts  of  smoothing  were  applied  to  the  wind 
field.  Data-void  areas  in  the  wind  field  were  filled  by  extrapolation. 
Three-point  and  nine-point  smoothing  operators  were  employed  to  smooth 
the  wind  field  prior  to  retrieving  the  perturbation  pressure  field. 
Smoothing  the  wind  field  produced  a  nonlinear  variation  of  the  central 
pressure  value  in  a  meso-low  with  respect  to  the  smoothing  parameter. 

The  central  pressure  value  decreased  as  the  smoothing  parameter  in¬ 
creased.  Pressure  centers  less  than  0.1  kPa  (1  mb)  in  absolute  magni¬ 
tude  tended  to  amplify  and  migrate  slightly.  False  pressure  centers 
created-  by  the  data-void  filling  technique  were  removed  when  a  smoothing 
parameter  greater  than  or  equal  to  0.4  was  employed.  Smoothing  para¬ 
meters  greater  than  or  equal  to  0.5  created  a  false  pressure  center. 
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1.  INTRODUCTION 


The  operational  meteorologist  is  given  the  responsibility  to  inform 
the  general  public  of  imminent  severe  weather  (e.g.,  tornadoes,  flash 
floods,  severe  thunderstorms,  damaging  winds,  damaging  hail).  Since  its 
inception,  the  conventional  weather  radar  has  been  the  meteorologist's 
primary  tool  for  detecting  and  monitoring  severe  weather  phenomena. 
However,  it  has  certain  limitations.  The  conventional  weather  radar  is 
only  capable  of  representing  the  spatial  distribution,  continuity  of 
movement,  and  intensity  of  the  precipitation  echoes  relative  to  severe 
weather  phenomena.  It  is  incapable  of  depicting  storm  kinematics, 
dynamics,  or  thermodynamics.  The  availability  of  some  of  this  informa¬ 
tion  might  improve  the  timely  issuance  of  weather  warnings  to  the  pub¬ 
lic,  decrease  false  alarms  by  improving  the  accuracy  of  warnings,  or 
provide  insight  into  the  nature  and  behavior  of  phenomena  peculiar  to 
convective  storm  activity.  Part  of  this  information  is  now  available 
from  Doppler  radars. 

With  the  advent  of  the  Doppler  radar,  meteorologists  are  now  able 
to  examine  the  internal  kinematics  of  convective  storms.  Single  and 
multi  pi  e-Doppler  radars  provide  velocity  data  from  which  the  internal 
kinematic  structure  is  readily  obtained.  In  recent  years,  it  has  been 
found  feasible  to  obtain  information  on  the  thermodynamics  and  dynamics 
of  convective  phenomena  through  the  use  of  the  four-dimensional  velocity 
information  obtained  from  the  Doppler  radar  by  solving  the  relevant 
dynamical  equations. 

The  citations  on  the  following  pages  follows  the  style  of  the 
Journal  of  the  Atmospheric  Sciences. 
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In  seeking  dynamic  and  thermodynamic  information  from  Doppler  data, 
particular  attention  has  been  given  to  the  tornadic  thunderstorm  and 
significant  insight  has  been  obtained  on  their  dynamic  structure.  Since 
the  method  of  retrieving  dynamical  information  from  kinematic  data 
promises  to  be  a  powerful  tool  in  our  analysis  of  convective  storms,  we 
need  to  carefully  analyze  the  implication  of  the  numerical  methods  used 
in  such  retrieval.  One  serious  problem  appears  to  be  the  noise  in  what¬ 
ever  wind  data  are  available.  This  problem  is  compounded  by  the 
presence  of  data-void  regions  due  to  the  absence  of  adequate  precipita¬ 
tion,  which  need  to  be  filled  by  some  extrapolation  or  interpolation 
procedure  which  itself  can  be  a  source  of  noise.  In  order  to  circumvent 
the  noise  problem,  some  kind  of  smoothing  scheme  is  applied  to  the  wind 
data  before  they  are  used  in  the  retrieval  scheme. 

This  study  will  examine  the  role  of  smoothing  on  the  retrieved 
perturbation  pressure  field  obtained  from  the  dual-Doppler  radar  data. 
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2.  BACKGROUND 

In  the  late  1970's,  as  indicated  above,  the  retrieval  of  thermo¬ 
dynamic  information  from  Doppler  radar  was  demonstrated  to  be  feasible. 
Initial  research  by  Bonesteele  and  Lin  (1978),  Hane  and  Scott  (1978), 
and  Gal-Chen  (1978)  showed  that  the  kinematic  information  from  Doppler 
radar  could  be  used  in  the  dynamical  equations  to  deduce  the  pressure 
distribution  associated  with  the  wind  field. 

Bonestelle  and  Lin  (1978)  worked  with  actual  storm  data.  They  com¬ 
bined  two  different  sets  of  single-Doppler  radar  data  from  two  tornadic 
thunderstorms  having  similar  characteristics  to  build  one  composite 
storm.  The  temperature  distribution  was  specified  around  the  boundary 
of  the  computational  domain,  and  the  (horizontal)  perturbation  pressure 
field  was  then  retrieved  from  the  composite  data  set.  In  their  retrie¬ 
val  scheme,  they  assumed  the  storm  to  be  in  the  steady  state.  Their 
results  indicated  that  perturbation  pressure  forces  were  important  for 
accelerating  moist  flow  toward  the  updraft  and  then  upward  into  it, 
while  these  forces  decelerated  the  downdraft  in  the  area  of  the  meso- 
high  and  induced  updrafts  along  the  gust  front. 

Hane  and  Scott  (1978)  and  Gal-Chen  (1978)  employed  synthetic  data 
generated  from  convective  cloud  models  to  test  their  retrieval  schemes. 
Gal-Chen  (1978)  introduced  pseudo-random  noise  into  the  test  data  to 
determine  how  it  affected  the  retrieved  thermodynamic  fields  (perturba¬ 
tion  pressure  and  temperature).  In  both  cases,  the  steady  state  assump¬ 
tion  was  examined  and  found  to  be  partially  invalid. 

Research  conducted  during  the  early  1980's  (Clark  et  al_. ,  1980; 

Hane  et  al . ,  1980)  continued  the  focus  on  the  error  caused  when 
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different  levels  of  noise  are  introduced  into  the  synthetic  data.  Time 
evolution,  advection,  sampling  time,  and  pseudo-random  errors  (noise) 
were  evaluated  by  simulating  radar  data  collection  through  convective 
cloud  models.  The  thermodynamic  retrieval  schemes  were  found  to  be 
sensitive  to  the  sampling  time  (frequency  of  observations)  and  white 
noise.  These  analyses  point  up  the  question,  already  asked  by  Hane  and 
Scott  (1978),  namely,  whether  a  filtering  of  Doppler-radar  data  suffi¬ 
cient  to  suppress  unrepresentative  noise  will  still  allow  the  calcula¬ 
tion  of  meaningful  results.  Another  problem,  related  to  using  observed 
data  with  retrieval  schemes,  was  identified  by  Hane  and  Scott  (1978). 

It  involved  the  filling  of  data  voids.  The  computational  schemes  re¬ 
quire  velocity  information  at  every  grid  point  in  a  three-dimensional 
Cartesian  domain,  whereas  observed  data  will  only  be  available  within 
irregular  boundaries.  They  also  noted  that  care  should  be  exercised 
when  interfacing  environmental  data  with  Doppler  data. 

Gal -Chen  (1978)  experimented  with  a  linear  smoothing  operator  to 
filter  noise.  From  Doppler  radar  experiments,  conducted  in  the  plane¬ 
tary  boundary  layer  to  determine  heat  and  momentum  fluxes.  Gal -Chen  and 
Kropfli  (1984)  found  that  a  "more  accurate  pressure  field"  resulted  if 
the  velocity  field  is  smoothed  prior  to  the  retrieval  of  the  pressure 
and  temperature  perturbations.  Brandes  (1984b)  also  applied  a  smoothing 
scheme  to  the  Doppler  radar-derived  velocity  field  prior  to  retrieving 
pressure  and  buoyancy  perturbations. 

With  it  being  apparent  that  some  form  of  smoothing  is  necessary. 


the  effort  of  this  study  will  be  to  examine  the  behavior  of  the  pertur¬ 
bation  pressure  field  when  different  amounts  of  smoothing  are  applied  to 
the  Doppler  radar-derived  wind  field. 


3.  THEORY 


As  mentioned  above,  retrieval  of  thermodynamic  information  from 
Doppler  radars  was  realized  in  the  late  1970's.  The  retrieval  procedure 
i.n  this  study  is  based  on  the  horizontal  momentum  equations  and  the 
anelastic  continuity  equation  (Pasken  and  Lin,  1982).  Because  it  is 
convenient  to  use  a  rectangular  computational  domain,  data  voids  adja¬ 
cent  to  the  observed  data  are  filled  via  an  objective  analysis  techni¬ 
que.  Since  the  data  generated  to  fill  in  the  data-void  areas  are  ficti¬ 
tious,  a  smoothing  scheme  is  necessary  to  decrease  excessive  gradients 
generated  in  the  vicinity  of  such  data  (Brandes,  1984a). 

a.  Perturbation  pressure  retrieval  scheme 

For  retrieval  of  the  perturbation  pressure,  the  assumption  of  a 
hydrostatic,  adiabatic  base-state  atmosphere  is  easily  justified.  In 
addition,  the  governing  equations  can  be  simplified  through  the  anelas¬ 
tic  approximation  (Ogura  and  Phillips,  1962).  Two  additional  assump¬ 
tions,  that  of  a  steady  state  and  a  frictionless  atmosphere,  while 
generally  invalid,  are  imposed  for  simplicity.  It  may  be  noted  that  the 
study  being  reported  has  rather  limited  objectives,  and  these  assump¬ 
tions  do  not  detract  from  these  objectives. 

With  the  above  assumptions,  the  governing  (quasi-horizontal)  equa¬ 


tions  are  as  follows: 


V 


V*  (PgC)  =  0  , 


P  =  Pg  +  P'  . 


where  C  is  the  three-dimensional  velocity  vector,  u  and  v  are  the  x  and 
y  components  of  the  horizontal  velocity,  respectively,  p  is  the  total 
pressure,  p^  is  the  environmental  hydrostatic  pressure,  p'  is  the  per¬ 
turbation  pressure,  is  the  constant  Coriolis  parameter,  v  is  the 
three-dimensional  gradient  operator,  and  is  the  environmental  density 
which  is  a  function  of  height  only. 

Following  the  method  of  Pasken  and  Lin  (1982),  (1)  and  (2)  are  put 
in  the  flux  form  as  follows: 
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By  taking  the  first  derivatives  with  respect  to  x  and  y  of  (5)  and 
(6),  respectively,  and  then  adding  them  together,  the  following  Poisson 
equation  for  perturbation  pressure  is  obtained: 


9  V  9  U\  ( 
9x  ■  9y'  "I 


9^(pgU^)  9^(p  UV)  9^(p  V^) 

- 2 - ^  2  3X3V  - V- 

9X^  9y‘^ 


(Cont. ) 


where  is  the  (horizontal)  two-dimensional  Laplacian  operator.  Since 
the  domain  used  is  too  small  to  be  free  from  the  effect  of  the  storm, 
the  lateral  boundary  conditions  imposed  are  of  the  Neumann  type,  namely: 


PgC'Vu  +  PgfpV  , 


=  -  p^C-Vv  -  . 


Since  u  and  v  wind  components  are  available  from  Doppler  radar,  the 
vertical  mass  flux  per  unit  area  can  be  obtained  from  (3)  with  no-slip 
surface  condition.  Perturbation  pressure  can  then  be  obtained  from  (7) 
with  boundary  conditions  (8)  and  (9)  since  the  forcing  terms  in  (7)  are 
known  functions  of  position.  It  may  be  noted  that  the  Coriolis  terms 
could  be  scaled  out  of  the  momentum  equations,  but  they  are  retained 
since  the  resulting  increase  in  complexity  is  small.  It  is  important  to 
note  that  three  dimensionality  enters  into  these  equations  through  the 
vertical  advection  terms. 


b.  Objective  analysis  for  filling  of  data-void  regions 


A  weighted  mean  recursive  interpolation  scheme  was  used  to  estimate 
grid  point  values  of  u  and  v  wind  components  in  the  data-void  regions 
(Cressman,  1960;  and  Trares  and  Das,  1982).  This  method  of  objective 
analysis  represents  a  quantity  at  a  grid  point  as  the  sum  of  the  weight¬ 
ed  values  of  data,  within  a  specified  radius  of  influence,  surrounding 
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the  quantity  to  be  interpolated.  With  this  procedure,  values  are  extra¬ 
polated  and  interpolated  into  the  data-void  areas  by  two  successive 
applications  of  this  scheme.  A  10  km  radius  of  influence  was  used.  The 
weight  function  is  dependent  on  the  influence  radius  and  the  distance 
between  grid  points  (resolution). 

The  interpolated  quantity  at  a  grid  point  is  defined  as  follows: 

_  M  M 

q(x,y)  =  Y.  5(d.>")q(x.y)  /  Y  c(d.r)  .  (10) 

m=1  m=l 

where  x  and  y  are  the  Cartesian  coordinates  of  a  grid  point,  r  is  the 
radius  of  influence,  d  is  the  distance  between  the  interpolated  quantity 
(^)  and  the  weighted  quantity  (q),  M  is  the  number  of  points  within  the 
radius  of  influence,  and  5  is  the  weight  function.  The  particular 
weight  function  employed  is  that  devised  by  Barnes  (1964).  It  is  de¬ 
fined  as  follows: 

C(d,r)  =  exp(-4d^/r^)  ,  (11) 

where  4  represents  a  reliability  factor  related  to  the  ability  of  the 
weight  function  to  represent  the  influence  of  the  data  within  the  radius 
of  influence.  As  with  any  weighted  mean  interpolation  scheme,  a  certain 
degree  of  smoothing  is  present  (Barnes,  1964;  Trares  and  Das,  1982). 

c.  Smoothing  scheme 

Three-point  and  nine-point  smoothing  operators  were  employed  to  ex¬ 
amine  how  different  amounts  of  smoothing  of  the  wind  field  affected  the 
perturbation  pressure  field.  The  (one-dimensional)  three-point  operator 


IS  defined  as  follows: 


*  f  (G,^,  *G,.,) 


(12) 


Application  of  this  operator  in  two  dimensions  (or  directions)  gives  the 
nine-point  smoothing  operator 


^‘^i+l.j+1  ^--l.j+l  ‘^i+l.j-l  ^i-l.j-1^/^ 


RG.^.  , 


(13) 


where  i  and  j  represent  coordinates  of  a  grid  point,  G.  .  is  the  origi- 
nal  value  of  the  quantity  to  be  smoothed,  G’l'^.  represents  the  smoothed 

»  »J 

quantity,  ij  superscript  means  smoothing  in  the  x-direction  then  in  the 
y-direction  with  the  (one-dimensional)  three-point  operator  (12),  a  is 
the  smoothing  parameter,  and  R  is  the  response  function. 

If  the  spatial  distribution  of  a  variable  can  be  represented  as  a 
harmonic  function,  then  the  response  function  obtained  by  inserting  a 
simple  harmonic  function  of  the  form 

G(x,y)  =  Aexp(i(kx  +  Hy))  ,  (14) 

into  (13)  gives  a  response  function  of 


R  =  (1  -  2asin^(^))  *  (1  -  2asin2(^))  ,  (15) 


where  A  represents  the  wave  amplitude,  Ax  and  Ay  are  the  grid  distances 


in  the  x  and  y  directions,  respectively,  and  k  and  £  are  the  wave  num¬ 
bers  in  the  x  and  y  directions,  respectively.  For  smoothing  parameters 
less  than  or  equal  to  0.5,  the  smoothing  schemes  do  not  affect  the  wave 
number  or  the  phase  of  the  original  wave,  and  are  independent  of  the 
first  direction  smoothed.  With  a  smoothing  parameter  greater  than  0.5, 
the  one-dimensional  smoothing  operator  will  produce  a  phase  shift  and 

the  two-dimensional  smoothing  operator  will  allow  shortwave  noise  to 

appear.  Fig.  1  is  a  graph  of  the  response  function  (R)  for  different 

values  of  the  smoothing  parameter  with  k  =  2"^/^  and  Ax  =  Ay  =  d.  For 

a  =  0.5  and  x  =  lid,  the  response  function  is  R  =  0.85.  This  means  that 
a  wave  with  a  wavelength  of  lid  will  lose  only  15%  of  its  amplitude  when 


smoothed.  On  the  other  hand,  for  the  same  smoothing  parameter,  a  wave 
with  a  wavelength  of  x  =  2d  will  be  completely  removed  (Shuman,  1958; 
Haltiner  and  Williams,  1980). 


4.  NUMERICAL  METHODS 


a.  Numerical  grid 

The  numerical  work  is  performed  on  a  standard  three-dimensional 
Cartesian  grid.  The  points  are  spaced  evenly  every  1000  m  in  the  hori¬ 
zontal  and  vertical  directions.  Only  two  vertical  levels  above  the  sur¬ 
face  are  considered  since  the  study  is  more  for  examining  the  soundness 
of  the  methodology  than  for  unravelling  storm  characteristics. 

In  the  difference  equations  presented  below,  the  grid  intervals  are 
Ax  =  Ay  =  Az  =  d  (=  1000  m).  The  customary  indices  i  and  j  are  used  for 
the  grid  points  in  the  x-  and  y-directions,  respectively,  and  the  sub¬ 
script  n  is  used  to  indicate  vertical  levels. 

b.  Finite-difference  equations 

(i)  Anelastic  Continuity  Equation.  The  vertical  mass  flux  per 
unit  area  was  obtained  by  vertical  integration  of  the  anelastic  continu¬ 
ity  equation  (3)  in  finite-difference  form  as  follows: 


‘  ‘^e'^i+1  ,j  ,n+l  '  ,n+l 


*  ^  '  Vi  J-l  ,n+l  >^2cl)/2)*d  (16) 


The  vertical  mass  flux  was  computed  at  the  first  two  vertical  levels 
above  the  bottom  boundary  of  the  computational  domain.  It  may  be  noted 
that  a  trapezoidal  form  is  used  in  the  vertical  with  a  centered  differ¬ 
ence  in  the  horizontal.  Vertical  velocity  is  treated  as  zero  at  the 
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lower  boundary. 

(ii)  Perturbation  Pressure.  Equation  (9)  is  put  in  finite- 
difference  form  as  follows: 


^i+l.j  ^  ^i-l.j  ^i,j+l  ^i,j-l  "  ^*^i,j 


'((u; 


2*u?  i  +  u?  ,  .)/d^ 


e  "“1+1, j  “  “1,j  “i-l,j 


-  (<“''e*')l+1,n+l  -  ('‘V’i+l.n-l  '  (“'>e“)i-l  .n+l  *  (%")i.l  ,n-l>''‘'‘ 

-  ((''V’j+l.n+l  -  (%“>j+l,n-l  -  (''V>j-l.n+l  *  <''»e“>J-l  .n-1 

where  P  represents  the  perturbation  pressure  in  the  above  finite- 
difference  formula. 


(17) 


(iii)  Poisson  Solver.  The  alternating  direction  implicit  method 
is  employed  to  solve  (17)  which  is  put  into  an  implicit-iteration  scheme 
as  follows: 


"’I.J  *  “‘‘"l+l-j  ■  '*"1.3  * ''l-l,J> 


t+1 


"“*<"1,3+1  -'*"1.3*  "1.3-1 3*8  '  ('8) 


where  B  is  an  acceleration  parameter,  t  is  an  iteration  level,  and  H-  . 

•  »J 

(forcing  function)  represents  the  right-hand  side  of  (7).  Equation  (18) 
is  solved  for  the  rows  of  the  numerical  grid  and  (19)  is  solved  for  the 
columns  of  the  numerical  grid  using  the  solution  obtained  from  (18), 
Equation  (18)  is  again  solved  with  information  obtained  from  (19).  This 
numerical  procedure  is  repeated  until  a  desired  level  of  tolerance  is 
obtained.  All  terms  in  this  scheme  and  in  the  forcing  function  are  com¬ 
puted  via  standard  centered  differences  (Gerald  and  Wheatly,  1984), 
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5.  DATA  AND  PROCEDURE 

a.  Data 

Dual-Doppler  radar  data  used  in  this  study  were  collected  by  the 
National  Severe  Storms  Laboratory  (NSSL)  on  2  May  during  project  SESAME 
(Severe  Environmental  Storms  And  Mesoscale  Experiment)  conducted  in 
1979.  NSSL  and  its  corresponding  network  of  Doppler  radars  are  located 
in  Oklahoma. 

On  the  morning  of  2  May  1979,  synoptic  charts  indicated  the  proba¬ 
bility  of  further  thunderstorm  development  from  pre-existing  nocturnal 
thunderstorm  activity  in  Oklahoma.  The  stability  indices,  vertical  wind 
shear,  frontal  position,  and  approaching  50  kPa  short  wave  trough  were 
conducive  for  possible  severe  thunderstorm  development  during  the  next 
24  h.  By  mid-afternoon,  isolated  thunderstorm  activity  developed  about 
175  km  northwest  of  Norman,  Oklahoma.  The  thunderstorm  activity  inten¬ 
sified,  and  two  large  cells  spawned  tornadoes.  Doppler-radar  data 
collection  began  around  1630  LST.  The  data  used  in  this  study  were 
collected  at  1651  LST  ( A1  berty  et  al_. ,  1979). 

The  basic  data  used  in  this  study  are  the  u  and  v  wind  components 
and  density  taken  from  a  vertical  sounding.  One  section  of  a  line  of 
tornadic  thunderstorms  associated  with  a  meso-circulation  was  examined 
in  this  study.  The  wind  component  data  have  already  been  assigned  to  a 
Cartesian  grid  spaced  at  1000  m  intervals  in  the  horizontal  and  vertical 
directions.  The  sounding  data  have  already  been  assigned  to  1000  m  in¬ 
tervals  in  the  vertical  direction.  Fig.  2  is  a  display  of  the  gridded 
wind  field,  derived  from  dual -Doppler  radars,  which  includes  general 
data-void  areas. 
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Fig.  2.  Gridded  wind  field  derived  from  dual-Doppler  radars, 
wind  vector  0.45  cm  in  length  represents  a  velocity  of  50  ms"'. 
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As  can  be  seen  in  Fig.  2,  the  samples  of  velocity  data  from  dual- 
Doppler  radars  have  some  inherent  problems,  one  of  which  is  that  infor¬ 
mation  is  only  available  within  irregular  boundaries.  In  addition,  the 
velocity  measurements  are  contaminated  by  a  considerable  amount  of  noise 
with  typical  errors  of  1-3  ms~^  or  more  in  areas  of  low  signal  intensi¬ 
ty.  Other  possible  problems  exist  with  the  data  used  in  this  study, 
namely,  those  due  to  the  location  of  storms  from  the  radars  (40-115  km), 
which  impact  on  the  resolution  of  low-level  phenomena  and  small-scale 
features  (Brandes,  1984b).  Fortunately,  these  do  not  seriously  affect 
the  present  study. 

The  nature  of  the  storms  and  data  collection  presents  wind  informa¬ 
tion  contained  within  irregular  boundaries.  To  avoid  an  irregular  com¬ 
putation  domain  and  for  convenience,  data-void  areas  in  the  observed 
wind  field  were  filled  prior  to  the  retrieval  of  the  perturbation  pres¬ 
sure.  The  data-void  filling  procedure  was  applied  to  minimize  forcing 
(i.e.,  excessive  gradients)  between  areas  with  data  and  areas  without 
data  so  that  a  rectangular  computational  domain  could  be  used.  All  data 
voids  in  the  wind  field  were  filled  by  extrapolating  and  interpolating 
values  into  the  data-void  regions. 


b.  Smoothing 


The  theory  and  properties  of  the  smoothing  operator  were  discussed 
in  Chapter  3.  The  three-point  and  nine-point  smoothing  operators  are 
utilized  to  smooth  the  Doppler  derived  wind  field  prior  to  retrieving 
the  perturbation  pressure.  This  was  done  to  suppress  noise,  to  suppress 
discontinuities  generated  by  the  data-void  filling  procedure,  and  pri¬ 
marily  to  study  the  response  of  the  perturbation  pressure  field  to  the 


smoothing  of  the  wind  field. 

The  boundary  winds  are  smoothed  with  the  three-point  smoothing 
operator,  and  the  rest  of  the  wind  field  is  smoothed  with  the  nine-point 
smoothing  operator. 

Different  amounts  of  smoothing  were  applied  to  the  wind  field, 
after  the  filling  of  the  data  voids  and  prior  to  retrieving  the  pertur¬ 
bation  pressure,  to  determine  how  the  smoothing  influenced  the  perturba¬ 
tion  pressure  field.  Values  of  the  smoothing  parameter  used  ranged  from 
0.0  to  1.0  with  0.1  increments. 

c.  Perturbation  pressure  retrieval 

The  scheme  for  the  retrieval  of  perturbation  pressure  from  dual- 
Doppler  radar  derived  wind  fields  has  been  described  earlier.  With  the 
u  and  V  wind  components  available  from  Doppler  radar  and  density  avail¬ 
able  from  a  vertical  sounding,  the  vertical  mass  flux  per  unit  area  was 
obtained  from  the  anelastic  continuity  equation.  This  makes  the  right- 
hand  side  of  (7)  a  known  function  of  position  so  that  the  Poisson  equa¬ 
tion  for  perturbation  pressure  can  be  solved. 

The  particular  numerical  technique  used  is  the  alternating  direc¬ 
tion  implicit  (ADI)  procedure.  In  implementing  this  numerical  proce¬ 
dure,  (18)  is  solved  by  sweeping  only  the  rows  of  the  numerical  grid  and 
(19)  by  sweeping  the  columns  of  the  numerical  grid,  using  information 
from  (18),  to  obtain  a  solution.  Equation  (18)  is  again  solved  with  in¬ 
formation  obtained  from  (19).  This  numerical  procedure  is  repeated 
until  the  change  in  the  solution,  between  successive  iteration  stages, 

has  reached  (or  become  less  than)  an  assigned  tolerance  which  was  set  at 
-5 

10  .  A  6  =  4.0  was  found  to  give  the  fastest  convergence  for 


streamfunction  retrieval. 


-5 

It  was  found  that  a  convergence  criterion  smaller  than  10  will 

cause  the  number  of  iterations  to  go  from  146  to  over  1000  using  the 

ADI  scheme  when  implemented  on  the  Harris  500  computer  in  the  single 

-5 

precision  mode.  For  the  problem  in  hand,  the  tolerance  of  10  was  con¬ 
sidered  adequate. 


6.  RESULTS 


The  perturbation  pressure  patterns,  associated  with  different 
amounts  of  smoothing  applied  to  the  wind  field,  have  been  contoured  and 
displayed  in  Figs.  3-13. 

Fig.  3  is  the  perturbation  pressure  pattern  associated  with  the  un¬ 
smoothed  wind  field.  The  meso-low  (-1.55  kPa)  pressure  area,  located  at 
X  =  8  km  and  y  =  11  km,  coincides  with  the  meso-circulation  seen  in  Fig. 
2;  this  is  in  agreement  with  previous  research  (Bonesteele  and  Lin, 
1978).  The  two  meso-high  pressure  centers  (0.549  and  0.554  kPa)  in  Fig. 
3  appear  to  be  due  to  extraneous  forcing  caused  by  the  data-void  filling 
procedure.  There  is  no  evidence  in  the  wind  pattern  that  supports  these 
perturbation  high-pressure  centers.  The  broad  low  pressure  area  to  the 
southeast  of  the  meso-low  pressure  center  is  characteristic  of  the  per¬ 
turbation  pressure  pattern  associated  with  meso-circulations;  this  is 
also  in  agreement  with  previous  research  (Bonesteele  and  Lin,  1978). 

Figs.  4-13  show  how  the  perturbation  pressure  field  responds  to  in¬ 
creased  smoothing  (increase  in  the  value  of  the  smoothing  parameter)  of 
the  wind  field.  The  false  perturbation-high  pressure  centers  are 
smoothed  out  except  for  one  located  at  x  =  21  km  and  y  =  8  km.  A  false 
pressure  center  appears  after  the  wind  field  is  smoothed  with  the  use  of 
a  smoothing  parameter  greater  than  or  equal  to  0.5.  This  is  evident  in 
Figs.  8-13  at  x  =  23  km  and  y  =  19  km.  This  pressure  center  continues 
to  amplify  as  the  wind  field  is  increasingly  smoothed  (Figs.  9-13).  The 
two  pressure  centers  located  at  x  =  2  km  and  y  =  8  km  and  x  =  1  km  and 
y  =  4  km  are  amplifying  and  decaying,  respectively,  in  amplitude  as  more 
smoothing  is  applied  to  the  wind  field. 


Fig.  3.  Perturbation  pressure  field  for  an  unsmoothed  wind  field. 
Units  are  0.1  kPa.  Contour  interval  is  0.2  kPa.  Relative  high  and  low 
pressure  centers  are  denoted  by  H  and  L,  respectively. 


Fig.  4.  Perturbation  pressure  field  for 
Contour  interval  is  0.2  kPa. 


■ig.  8.  Perturbation  pressure  field  for  a  =  0.5.  Units  are  0.1 
Contour  interval  is  0.2  kPa. 


Fig.  13.  Perturbation  pressure  field  for  a  =  1.0.  Units  are  0.1 
Contour  interval  is  0.2  kPa. 
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From  examining  Figs.  4-13,  it  is  apparent  that  a  smoothing  parame¬ 
ter  greater  than  or  equal  to  0.5  will  create  a  false  pressure  center  and 
cause  the  smooth  isobar  pattern  obtained  from  lighter  smoothing  to  be¬ 
come  noisy  again.  Also,  values  less  than  0.1  kPa  (absolute  value)  tend 
to  amplify  and  migrate  slightly.  This  is  apparent  from  the  pressure 
value  located  initially  at  x  =  3  km  and  y  =  8  km  in  Fig.  3.  In  Fig.  13, 
this  pressure  center  has  moved  to  x  =  4  km  and  y  =  8  km. 

Fig.  14  is  a  graph  of  the  center  pressure  value  of  the  meso-low  for 
different  values  of  the  smoothing  parameter,  normalized  by  the  central 
pressure  value  obtained  from  the  unsmoothed  wind  field,  plotted  against 
the  smoothing  parameter.  The  relationship  is  nonlinear  but  appears 
close  to  being  linear.  This  may  be  due  to  the  fact  that  the  smoothing 
operator  is  linear  whereas  the  forcing  for  the  perturbation  pressure  in 
(7)  is  nonlinear. 

Figs.  15  and  16  are  the  divergence  and  vorticity  fields  associated 

with  the  unsmoothed  wind  field,  respectively.  Values  of  divergence  and 

-2  -1 

vorticity  are  in  units  of  10  s  .  An  approximate  theoretical  estimate 
of  vorticity  in  the  meso-low,  assumed  to  be  a  Rankine  vortex,  was  made 
to  be  consistent  with  the  pressure  value  in  the  meso-low  for  the  un¬ 
smoothed  wind  (Fig.  3).  The  value  of  the  vorticity  obtained,  assuming  a 

-2  -1 

radius  of  3000  m,  is  2.6  x  10  s  .  In  Fig.  16,  the  value  is  3.9  x 
-2  -1 

10  s  .  The  theoretical  value  for  vorticity  agrees  fairly  well  with 
the  actual  value  of  vorticity.  The  difference  may  be  attributed  to  the 
divergence  field  associated  with  the  meso-circulation  (Fig.  2  and  Fig. 
15)  since  the  wind  field  of  the  Rankine  vortex  is  non-divergent. 

The  0.644  kPa  meso-high  pressure  in  Fig.  3,  x  =  22  km  and  y  =  8  km, 
has  some  support  from  the  divergence  field  displayed  in  Fig.  15. 


SMOOTHING  PARAMETER  (a) 


Fig.  14.  Effect  of  smoothing  on  retrieved  perturbation  pressure  at 
the  center  of  the  principal  meso-low.  is  perturbation  pressure  ob¬ 
tained  with  the  smoothing  parameter  (c).  P^  is  perturbation  pressure 
obtained  from  unsmoothed  wind  field. 


Fig.  16.  Vorticity  field  associated  with  the  upsmopthed  wind 
field.  Units  are  10"^  s"' .  Contour  interval  is  10"^  s" ' .  Positive  and 
negative  vorticity  centers  are  denoted  by  H  and  L,  respectively. 


However,  the  divergence  value  in  the  general  location  is  very  weak 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


a.  Conclusions 


Smoothing  the  wind  field  removes  most  false  pressure  centers 
generated  by  extraneous  forcing  from  the  data-void  filling  procedure. 
However,  values  of  the  smoothing  parameter  greater  than  or  equal  to  0.5 
created  a  false  pressure  center.  This  may  be  due  to  the  change  of  phase 
caused  by  the  three-point  smoothing  operator  used  to  smooth  the  boundary 
winds.  Increased  smoothing  of  the  wind  field  caused  pressure  centers  of 
less  than  0.1  kPa  (absolute  magnitude)  to  amplify  and  migrate  slightly. 

A  smoothing  parameter  of  0.4  will  remove  false  pressure  centers  genera¬ 
ted  by  the  data-void  filling  technique  without  introducing  more  short¬ 
wave  noise.  Increased  smoothing  with  a  smoothing  parameter  greater  than 
or  equal  to  0.5  caused  the  smooth  isobar  pattern  obtained  from  lighter 
smoothing  to  become  noisy  again. 

Based  on  the  perturbation  pressure  field,  the  wind  pattern  is  high¬ 
ly  ageostrophic  as  would  be  expected  by  the  scale  of  the  storm.  In  the 
meso-low,  the  wind  pattern  gives  the  suggestion  of  cyclostrophic 
balance.  However,  since  cyclostrophic  flow  is  nondivergent,  the  flow 
deviates  considerably  from  cyclostrophic  balance  because  the  wind  field 
in  the  meso-low  is  divergent. 


b.  Recommendations 


(i)  Steady  State.  If  possible,  the  steady  state  assumption  should 
be  discarded  and  the  local  time  derivatives  evaluated.  However,  this 
requires  short  radar  scan  periods  over  the  same  volume  sample  (Hane  and 
Scott,  1978;  Hane  et  al . ,  1980;  Clark  et  al . ,  1980;  Gal-Chen  and 


Kropfli,  1984;  Das,  1985). 


(ii)  Inclusion  of  Friction.  A  parameterization  of  friction  should 
be  included  in  the  equations  since  the  flow  is  not  inviscid,  and  kinetic 
energy  dissipation  needs  to  be  accounted  for  in  the  equations  of  motion 
(Das,  1985). 

(iii)  Three-Dimensional  Smoothing  and  Retrieval.  Because  the  wind 
field  associated  with  convective  Phenomena  is  a  three-dimensional 
entity,  the  smoothing  scheme  applied  to  the  wind  field  and  the  retrieval 
of  perturbation  pressure  needs  to  be  accomplished  through  three- 
dimensional  equations  (Das,  1985). 
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